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abstract 

\\ 

A concept  is  presented  which  facilitates  the  production  of  computer 
generated  plots  of  families  of  curves  in  either  2 or  3 dimensions. 

Specific  computer  examples  provided  are  intended  to  serve  as  master  programs 
from  which  a wide  variety  of  general  purpose  scientific  plotting  applications 
can  be  performed  with  a minimum  amount  of  effort. 
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SECTION  1. 

INTRODUCTION 

For  many  years  digital  computers  have  been  used  to  implement  math  models 
of  various  physical  systems.  However,  once  the  delight  of  obtaining  correct 
numerical  printouts  from  mammoth  digital  calculations  has  passed,  it  has  often 
been  an  additional  burden  to  effectively  and  inexpensively  display  the  results 
in  forms  such  that  the  scientist  can  gain  additional  insight  into  the  behavior 
of  the  system  he  is  studying.  In  order  to  alleviate  this  problem  there  have  been 
numerous  digital  plotting  systems  developed  which  are  capable  of  producing  any 
graph  imaginable  by  the  scientist.  However,  all  too  frequently  the  scientist 
is  hampered  by  various  roadblocks  which  usually  exist  in  such  endeavors.  First 
of  all,  the  lead  time  required  to  program  new  plot  configurations  can  often  be 
on  the  order  of  days.  This  lag  time  in  itself  can  be  sufficient  to  discourage 
the  innovative  scientist,  without  even  considering  the  manpower  costs  involved. 
This  paper  is  a description  of  work  performed  by  the  author  in  order  to  minimize 
these  bottlenecks  and  to  maximize  responsiveness  to  the  needs  of  scientific 
projects.  This  work  has  been  performed  by  the  NUC  Transducer*  and  Array 
Systems  Division;  however,  its  generality  should  be  applicable  to  a wide  variety 
of  other  scientific  projects. 

In  order  to  implement  the  concept  to  be  demonstrated  in  the  paper,  the 
following  two  items  are  required: 

1.  A highly  responsive  Computer  Center  with  plotting  hardware,  and 

2.  A high  level  user-oriented  graphics  computer  language. 

At  NUC  the  items  are  satisfied  by  the  following: 
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(a)  either  the  UNIVAC  1230/490  systems  with  their  associated  Calcomp  565 
and  1136  plotters,  or  the  UNIVAC  1108  and  its  associated  Calcomp  663,  and 

(b)  the  Fortran  software  system  DISSPLA*  as  implemented  on  the  NUC  UNIVAC 
computers  mentioned  above. 

The  concept  which  serves  as  the  basis  for  the  examples  to  be  shown  involves 
a data  storage  scheme  well  suited  to  the  problem  of  either  producing  families 
of  curves  in  2 dimensions,  or  producing  the  same  information  as  a perspective 
view  in  3 dimensions.  By  combining  the  simplicity  of  the  data  storage  scheme 
with  the  many  powerful  features  of  DISSPLA  it  has  been  possible  to  produce 
highly  versatile  plotting  samples  which  can  be  modified  for  various  applications 
with  a minimum  amount  of  effort. 

The  computer  examples  to  be  presented  in  this  paper  are  as  follows: 

1.  PLT  FAM  for  plotting  families  of  curves  in  2 dimensions 

2.  PLT  SRF  for  3 dimensional  plotting  of  either  families  of  curves  or 
surfaces 

3.  PLT  PTN,  a modified  version  of  PLT  FAM  specifically  oriented  towards 
polar  plot  applications. 

For  each  example,  the  data  storage  scheme  will  be  outlined,  the  program  op- 
tions will  be  explained,  and  the  associated  computer  printout  and  plots  will  be 
included. 


SECTION  2. 

COMPUTER  PROGRAMMING  CONSIDERATIONS 


The  programs  have  been  developed  on  the  NUC  UNIVAC  1230  computer  using  ASA 
Fortran  with  two  exceptions,  SFLD  and  ENCODE.  In  order  to  run  the  programs  on 
the  NUC  UNIVAC  1108,  several  Fortran  statements  (marked  in  the  decks  by  having 
"1230"  in  columns  73-76)  must  be  changed  as  shown  in  Fig.  2.1.  Similar  chanqes 
would  be  required  on  any  other  computer  system  using  the  DISSPLA  package. 

Actually  the  statements  requiring  SFLD  (FLD  on  the  1108)  can  be  eliminated, 
if  desired,  since  they  are  used  only  to  provide  a safety  factor  in  case  the 
user  fails  to  terminate  one  of  the  self-counting  title  strings  properly. 

Even  though  all  three  examples  were  originally  developed  on  the  1230, 
the  PLT  SRF  example  will  be  for  the  1108,  since  the  current  version  of  the 
DISSPLA  3D  software  is  too  large  for  the  1230  when  using  all  of  the  options 
demonstrated.  PLT  SRF  can  be  run  on  the  1230  if  several  of  the  options  are  deleted. 

In  order  to  implement  new  applications  the  user  must  first  of  all  select 
a data  input  medium.  Although  this  could  involve  direct  computation  of  the 
function  to  be  plotted,  it  is  frequently  advisable  or  necessary  to  load  pre- 
computed values  from  tape  or  cards  as  shown  in  the  examples.  Once  the  input 
medium  has  been  selected,  the  user  must  then  tailor  the  data  input  section  of 
the  CONTROL  ROUTINE  to  the  needs  of  his  problems  and  also  select  the  appropriate 
options  for  his  needs  from  one  of  the  three  examples. 


C***********  NUC  UNIVAC  1230  STATEMENTS  ************* 


UATA  [J  wRD  P C /16/  C WORDS  PER  CARD  IMAGE  ON  UNIVAC  1230 

DATA  ID0LAK/5H  $ / 1230 

CALL  SFLD(  24#  6#  I HEAD ( 12  # J) # IDOLAR  ) 1230 

200  FORMAT ( 16A5)  1230 

201  FORMAT ( IX  * 16A5)  1 1230 

CALL  SFLD ( 24 * 6#  IHD  LPT(6)»  IDOLAR)  1230 

CALL  SFLD(24»  6#  IHD  MEM(6)»  IDOLAR)  1230 

CALL  SFLD ( 24#  6#  IHDFAM ( 12  #N) # IDoLaR  ) 1230 

common  / cm  plot  / i bufi  200 ) set 

DATA  NBUF  / 200  / SET 

DATA  I PLT  UN  / 2 / SET 

DATA  ID0LAR/5H  S / 1230 

CALL  SFLD<  24#  6#  IHDFAM ( 12 » N ) » IDOLAR  ) 1230 

4990  FORMATS*  COMPLETED  PLOT  NO.*#  13#  ♦ ***  *»  12A5  / 1230 


C***********  fJUC  UNIVAC  1108  STATEMENTS  ************* 


DATA  N WRD  P C /14/  „„RDS  PER  CARO  IMAGE  On  UNIVAC  1108 

DATA  ID0LAR/6H  $/  1108 

FLD(30»  6#  I HEAD(12#J)  ) = I DOLAR  1108 

200  FORMAT<13A6#A2)  1108 

201  FORMAT ( IX  # 13A6.A2)  1108 

FLD(30»  6#  IHD  LPT(6)  ) = I DOLAR  1108 

FLD ( 30 » 6#  IHD  MEM ( 6 ) ) = I DOLAR  1108 

FLD (30 » 6#  IHD  FAM<12»N>  ) = I DOLAR  1108 

COMMON  / CM  PLOT  / I BUF(2000)  SET 

UATA  NBUF  / 2000/  SET 

UATA  I PLT  UN  / 8 / SET 

DATA  ID0LAR/6H  %/  1108 

FLD ( 30 » 6#  IHD  FAM(12»N)  ) = I DOLAR  1108 

4990  FORMAT ( / * COMPLETED  PLOT  NO.*#  13#  * ***  '»  12A6  / 1108 


FIG  2.1  - FORTRAN  STATEMENTS  TO  BE  CHANGED  FROM  123o  VERSION 

TO  PRODUCE  1108  VERSION  OF  PLT  FAM 
( SIMILAR  REQUIREMENTS  FOR  PLT  SRF  AND  PLT  PTN  ) 
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SECTION  3. 


PLT  FAM 


PLT  FAM  is  designed  for  plotting  families  of  curves  in  2 dimensions  with 
the  following  features: 

1.  Option  for  user  selected  axis  scales  or  automatic  self  scaling  via  DISSPLA 

2.  User  supplied  headings,  axis  labels,  and  legend  labels 

3.  Choice  of  linear,  log-log,  log-linear,  or  linear-log 

4.  Option  to  combine  two  consecutive  families  of  curves  onto  the  same  grid 

The  data  storage  scheme  for  the  X,  Y data  points  is  illustrated  in  Fig.  3.1, 
where  NFAMLY  is  the  number  of  families  of  curves,  MEMBER  is  the  number  of 
members  per  family,  and  LPTS  is  the  number  of  points  per  member  of  the  family. 
Comments  at  the  beginning  of  the  PLT  FAM  subroutine  explain  the  array  sizes  required. 

Detailed  line  by  line  descriptions  of  the  programs  will  not  be  given  since 
an  attempt  has  been  made  to  provide  sufficient  comments  within  the  Fortran 


programs . 


CERAMIC 


6,1971  [J 
LESS  HERD  0 "I 

a - 1 

IE  E0-5500+-T 
x -T 


I k \i  V M 

£ O CD 


LEGEND 
o-Tl  .0200 
a -T2  . 0292 

SEADUCER  RUN  80122  MAY  6,1971  + -T3  .0400 

IP  2.3  DESIGN  RUN  WITH  LOSSLESS  HEAD*  I™ 

TRIAL  DESIGN  GENERATING  SCHEME  P0-5500v -T6  !o800 

® -T7  .1000 

TAIL  LENGTH 


CERAMIC  AREA  (C2  THRU  C9) 


{ 


SECTION  3.2 

PLT  FAM  program  listing 
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JOB  60l  00.04  L MCCLEARY 

woUNT  M3S=PL0T 

fortpam 

***********  PLT  FAM  CONTROL  ROUTTNF  * * * 
COMMON  / / XDATA ( 30),  YDATA ( 30,12*6) 

COMMON  / C WMF  AO  / IHDLPT  (6)  » IHDMFM<6),  IHDFAM(1?,  6)  , 
COMMON  / NMMS  / LPTS,  MEMBER  , NFAMLY , NHLTNS 
Common  / CPL  LRL  / DTA  LBL(12»2),  IPK  LRLM2,12> 
DIMFNSION  I CHK(l6> 

OATA  fj  WRO  P C /1 6/  O WORDS  PPR  CARO  IMAGF 

DATA  I00LAR/5H  % / 

OATA  IN  UNIT  / S / 


MCCLFARY  PLT  FAM  R01220SBN 


IHFAO ( 12  * 3) 


IINIVAC 


1 230 
1230 


C**********  USFR  SFTS  NHLINS  = NUMRFR  OF  HEADING  LINES  TO  RF  LOAOFO 
NHL I MS  = 3 

C**********  LOAD  THFAO  POR  PLOT  HEADING^ 

DO  4100  J = 1» NHL INS 

RF AO ( IN  UNIT,  200  ) ( IHFAO ( I , J) , T = l»12) 

CALL  SFLD«  24,  6,  I HEAD(12,J),  IOOLAP  ) 


4100 


C********** 

READ  ( 


CALL  SFLD( 

PRINT 

CONTINUE 

*****  READ 


1230 


*****  READ  and  p 
READ ( IN  UNIT,  2 
PRINT  2 

READ ( IN  UNIT,  2 
PRINT  2 

FORVAT(  16A9) 
FORMAT (IX*  16AS) 


print 
200  ) 
201  , 
200  ) 
201  , 


( IHEAO ( I » J) * 1 = 1 » 12) 

iMY  EXTRA  INFORMATION 
(I  CHK(I) » 1 = 1 » N WRD 
(I  CHK(T) , 1 = 1 » N WRD 
(I  CHKCT),  T=1»N  WRD 
(I  CHK(T) » 1=1 »N  WRD 


IF  DESIRED 
P C) 

P C) 

P C) 

P C) 


1230 

1230 


C**********  USFP  LOADS  INPUT  DATA  FOR  LPTS*  MEMRFR,  NFAMLY, 


DTA  lpl 


POUNDS 
LPTS  = 

mempfr 

NFAMLY 
DO  9 M 
DO  P L 


= 2.20462 


1 »MEMRFp 
1 , LPTS 


PRINT  101* 

1 

ino  format< 
ioi  format ( ix , 

DO  10  N = 
IF  ( N .EQ 


REA0(5» 100)  DT ALnL ( M» 1 ) » DTALRL(M»2)»  XDATA (L) * 
(YDATA(L»M,N) , n=i»nfamly) 

NT  101*  DTALRL(M,1)»  DTALRL(M»2>»  XDATA (L) * 
(YDATA(L,M,N) , N=1»NFAMLY) 

MAT(  IX, A2,  4X,3F9.6*3X,FP.4,FR.5,3F10.4) 

MAT ( IX , IX, A2,  4X»3FR,6»3X»F°,4»Fft,S»3E10.4) 

10  N = 1, NFAMLY 

( N .EQ.  2 ) YDATA(L*M,N)  = Y DATA (L»M»N)  * PO 


DATA(L»M,N) 


POUNDS 


10  CONTINUE 

8 CONTINUE 

9 CONTINUE 


C**********  LOAD  I HD  LPT  FOR  X AXIS  TITLF 

READ ( IN  UNIT,  2oO  ) ( I HD  LPT ( I ) , 1=1,6) 

CALL  SFLO ( 24 , 6,  I HD  LPT(6>»  IDOLAR ) 

PRINT  20 l , ( IHn  LPT ( I ) , T=l,6) 

C**********  LOAD  I HD  MEM  FOP  LEGEND  TITLE 

RFAD ( IN  UNTT,  200  ) ( I HD  MEM(I),  1=1,6) 

CALL  SFLD124,  6,  THD  MEM(6) , IDOLAR) 

PRINT  201  » ( I HD  MEM(I),  T=1 ,6) 

C**********  LOAD  I HD  FAM  FOR  INDIVIDUAL  PLOTS  ( Y 
DO  4000  N = 1, NFaMLY 

REAO(  IN  UNIT,  200  ) ( IHD  FAM ( I , N) , 1=1,1?) 

CALL  SFLD<  24,  6,  IHDFAM ( 1 2 ,N) » IDOLAR  ) 
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AXIS  ) 


— 


PRINT  201  » (IHn  FAM<I»N) * T=l*l?) 

4000  continue: 

c * * * * * INPUT  IS  NOW  LOADED  ***** 

CALL  PLT  FA* 

C * * * * * terminate  fnwire  PLOT  RUN  ***** 

CALL  PLOT ( 0 , » 0.,  999) 

STOP 
END 


ooonooooono 


SUBROUTINE  PLT  PflM 

C * * * * * PLT  FAM  ( USING  DISSPLA  PACKAGE)  **************** 

c * * * * * general  program  to  plot  families  of  curves  with  legends  * * * * 

c*******  PROGRAMMED  PY  L.  E.  MCCLEARY  Nile  CODE  601  ******* 

*****  DIMENSION  AS  FOLLOWS  ***** 

*****  WHERF  L=LPTS»  MrMEMPFR#  NrNFAMLY,  HsMHLTWS 

* * * LOAD  COMMONS  / /»  /CMHFAD/ » /NUMS/»  AND  ARRAY  DTA  LRL  IN  CONTPOL 

* * * LOAD  COMMONS  /P  CMBIN/»  /PL  EORM/r  TN  THIS  OOUTINF  IF  NFEDEO 

COMMON  / / XDAT A ( L ) * YnATA (L»M»N) 

COMMON  / CmhEAD  / IHDLPT(6)»  THDMFM(6)»  INGRAM!!?,  N)  , IHFADC1?»H) 
COMMON  / MUMS  / LPTSr  MFMREP#  NF AMLY » MHLTNS 
COMMON  / CPL  LRL  / DTA  LBL ( M.2)»  I»K  LPL(  M,12) 

COMMON  / NO  LOAD  / YMIN(M)»  YmAX(N),  YPLOT(L) 

COMMON  / P CMRIN  / N CMBIN(N) 

COMMON  / PL  FORM  / N FORM(N)  0 LIN-0»  I.OG-1 1 XLG-2*  YL6-3 

COMMON  / / XDATA(  -AO).  YDATA(  30*12#6) 

COMMON  / CmweaD  / IHDLPT(6)»  IHDMFM(6>»  IHOFAM(l?t  6)  * IHF  AD  ( 1 2 » 3) 
COMMON  / MUMS  / LPTS*  MEMBER*  NFAMLY#  NHL INS 
COMMON  / CPL  LRL  / DTA  LRL(12»2),  IPK  LRL ( 1 2» 1 2) 

COMMON  / NO  LOAD  / YMIN(6)»  Y MAX<6>»  Y PLOT(30) 

COMMON  / P CMRTN  / N CMBIN<6> 

COMMON  / PL  FORM  / N F0RM(6)  D LIN-0#  LOG-1 » XLG-2#  YLG-3 

Common/  OQFYTR  / XOR # XSTFP,  YOR,  YSTEP # X AXIS#  Y AXIS#  LX#|.Y»TT 
DIMENSION  iw  LRL  M ( 4 ) 


COMMON  / CM  PLOT  / I RUF ( 200)  FrT 
DATA  NP.UF  / 200  / SFT 
DATA  IFTF/O/ 

DATA  I PLT  UN  / ? / SFT 
DATA  IDOLAR/SH  $ / 1230 


PRINT  1 
1 FORMAT ( 1 HI ) 

C * * * * * INITIALIZE  plot  puffer 

IFCIFTF  .EO.  0)  CALL  PLOTS ( IRUF»NRUF»  I PLT  UN) 

IFTFzl 

C CALL  F ACTOR ( . 85)  0 REDUCE  fl  1/2  RY  H TO  VIEW  GRAPH  ST?f 

C * * * * * FIND  max  MIN  OVER  ALL  MEMBERS  FOP  EACH  F AMTLY 
QO  GOO  N = 1#  nfamly 

Y MlN(N)  = 1.E30 

Y MAX (N)  = -1.E30 

C**********  USER  SETS  N FORM ( N)  =0»1#2»3  FOR  LlN>  l.OG»  XLG  OR  Yl  G PLOTS 
N FORM(N)  = 0 

C**********  USER  SETS  N CMRIN(N)  = 1 TO  PLOT  FAMLIES  N AND  N+l  ON  SAME  GRID 
N CMplN(N)  = 0 

500  CONTINUE 

DO  3000  L = 1 » LPTS 
DO  2000  M = 1»  MEMBER 
DO  1000  N = 1»NFAMLY 

IF  ( Y DATA (L»M»m)  .LT.  Y MIN(N)  ) Y MTN(N)  = Y DATA (L»M»N) 

IF  ( Y DATA(L»M»n)  .GT.  Y MAX(N)  ) Y MAX(N)  = Y PATA(L*M,M) 
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1000  CONTINUE 
2000  CONTINUE 
3000  CONTINUE 

C********  USE  nlSSPLA  PACKAGE  ******** 

C * * * * * LOOP  ON  M FPOM  1 TO  MFAMLY  ***** 

NTH  PLT  = 0 
N = 0 

4909  N = N + 1 

NTH  PLT  = NTH  PLT  ♦ 1 

IF  ( N CMBIN(N)  .FQ.  0 ) 00  TO  4101 

IF  ( Y MIN ( M+l ) .LT.  Y MJM(N)  ) Y MIN(*')  = Y MTN(N+1) 

IF  ( Y MAX (N+l ) .0 T.  Y MAX ( N)  ) Y MAX (M)  = Y MAX(N+1) 

DO  4100  N HEAD  = 1 , f> 

IHDFAM(NHEAD+6»N)  = IH0FAv(NHEAD»N+1) 

4100  CONTINUE 

CALL  SFLD(  24,  ft,  IHDFAM( 12,N) , TDOLAR  ) 1230 

4101  continue 

CALL  PGNPL ( -NTH  PLT  ) 

M = 1 

DO  4400  L = 1,LPTS 
4400  Y PLOT ( L ) = Y DATA(L,M,N) 

C**********  USEP  GETS  I SCALE  = 0 FOR  SFLF  SCALING 
C**********  USFR  SETS  T SCALF  = 1 FOR  USFR  SUPPLIED  SCALES 
I SCALE  = 1 
I SCALF  = 0 

IF  ( I SCALF  .FO.  0 ) 00  TO  4440 

c**********  for  user  Supplied  scales, 

c**********  USER  sets  I AXIS’ Y XXX S»  XOP,YOR,  XSTFP, YSTFP  FOR  f ACH  FAMILY 
GO  To  4445 
4440  CONTINUE 

C * * * * * SFLF  SCALING  TO  RETURN  VARIABLES  IN  DISSPLA  COMMON  /QQFXTP/ 

CALL  HLPPLT 

CALL  X RANGE ( X DATA ( 1 ) , X DATA (LPTS)  ) 

CALL  Y RANGF ( Y mTN(N),  Y MAX(M)  ) 


IF 

( 

N 

FORM { n ) 

,F0. 

0 

) 

CALL 

LINPLT ( 

xdata, 

YPLOT, 

LPTS 

) 

IF 

( 

N 

FORM(N) 

,FQ. 

1 

) 

CALL 

LOGPLT ( 

XDATA, 

YPLOT, 

LPTS 

) 

IF 

( 

N 

FORm(N) 

• EO  * 

2 

) 

CALL 

XLGPLT ( 

XDATA, 

YPLOT, 

LPTS 

) 

IF 

( 

N 

FORM (N) 

.eg. 

3 

) 

CALL 

YLGPLT ( 

XDATA, 

YPLOT, 

LPTS 

) 

4445  CONTINUE 

C * * * * * SCALES  ARE  NOW  AVAILABLE  FOR  PLOTTING 
ITITLF  = 1 

IF  ( X AXIS  ,GT.  Y AXIS  ) 7TITI.E  = -1 

CALL  TITLE ( I H ,ItTTLE*  IMDLPT»100*  THDFAMt 1 ,M) ,100,  X AXIS»Y  AXIS) 


IF 

( 

N 

FORM (N) 

» FO  . 

0 

) 

CALL 

GRAPH 

(XOR, 

XSTFP, 

YOR» 

YSTFP) 

IF 

( 

N 

FORM(N) 

• EG. 

1 

) 

CALL 

LOGLOGIXOP, 

XSTFP » 

YOR* 

YSTFP> 

IF 

( 

N 

FORM(N) 

.FO. 

2 

) 

CALL 

X LOG 

(XOR, 

XSTFP » 

YOR» 

YSTp°) 

IF 

( 

N 

FORM(N) 

.FO, 

3 

) 

CALL 

Y LOG 

(XOR, 

XSTFP » 

YORf 

YSTFp) 

DO  4500  J = 1 , NHL INS 

4500  CALL  HEAOINI  IHEaD(1#J)*  100,  NHL  IMS) 
CALL  FRAME 
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i. 


CALL  GPIO<  1 ) 


C CALL  SPLINE 

CALL  MARKER!  v ) 

I MARK  = 0 
I MARK  = 1 

CALL  CURVE ( X DATA . Y PLOT.  LPTS,  I MARK  1 

C**********  USER  SELECTS  LEGEND  FORMAT  FOR  FNCOOF 
4 FORMAT(A2*  R6.4.  *$') 

EMCODF ( 20 » 4.  TH  LRL  M)  DTALRL (M* ) ) » OTALPI  (M,?) 

CALL  L INES ( TH  LPL  M » IPK  LRL*  M ) 

IF  ( VEMREP  .EO.  1 ) GO  TO  49pl 
DO  4900  M = 2 » mEmREP 
DO  4600  L = l.LPTS 
4600  Y PLOT ( L ) = Y DATA(L.M.N) 

CALL  MARKER (M) 

CALL  CURVE ( X DATA.  Y PLOT.  LPTS.  I MARK  ) 

ENCODE  ( 20 » 4,  IH  LRL  M)  DTALRL  (M.  1 ) » PT  ALRl.  ( M»  2 ) 

CALL  LINES ( TH  LBL  M.  IPK  I RL.  M ) 

4900  CONTINUE 

4901  CONTTNUF 

C**********  USER  SETS  n LEGND  = 0 FOR  NO  LEGEND 
C**********  USER  SETS  N LEGND  = 1 FOR  LEGEND 
N LFGND  = 0 
M LFGND  = 1 

IP  ( N LFGND  ,FO.  0 ) GO  TO  4960 
XLEGrjD  = X AXIS  - 1. 

IF  ( X AXIS  ,LT.  6.  ) X LFGND  = X AXIS  + .25 
YLEGMD  = Y AXIS 

IF  ( Y AXIS  ,GT.  R.  ) Y LFGND  = Y AXIS  - ,R 

IF  ( ITITLF.NE.l  .AND.  YAXTS.GT.S.  > YLEGND  = YAXI?-.? 

CALL  LEGEND ( IPK  lBL*  MEMRFR,  XLEGND*  VLE«N"  ) 

CALL  MESSAG ( I HD  MEM.  100.  XLEGND.  YLFGND-.2  ) 

4960  CONTTNUF 

IF  ( N CMRIN(N)  .EG).  0 ) GO  TO  49R0 
DO  4970  Ms  1 » MFMBER 
DO  4965  L = 1.  LPTS 
4965  Y PLOT (L)  = Y DATA (L*M*N+t ) 

CALL  MARKEP ( M ) 

CALL  CURVE ( X DATA.  Y PLOT.  LPTS*  I MARK  ) 

4970  CONTINUE 
49R0  CONTINUE 

CALL  ENDPL ( -NTH  PLT  ) 

PRINT  4990.  NTH  PLT.  (IHO  FAM(T»N)*  1=1*12). 

P XOR.XSTEP.X  AXIS*  YOR * YSTFP. Y AXIS 

4990  FORMAT (/*  COMPLETED  PLOT  NO.*.  13*  ' ***  **  12A5  / 

F • X*.  2F10.4*F4.1*  » Y« * 2F10.4.F4.1  /> 

IF  ( N CMBIM(N)  .EQ.  0 ) GO  TO  5000 
N = N + 1 

5000  IF  ( N .LT.  NFAMlY  ) GO  TO  4909 
RETURN 
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I. 


SEND 

SDATA 

SEADUCEP  PUN  80122  MAY  6rl971$  LIKF  00089  Q=50 
IP  2.3  DESIGN  PUN  WITH  LOSSLESS  HEAPS 
TRIAL  DESIGN  GENERATING  SCHEME  F0=5500  % 


CHARGE 

TO  16000501  ARRAY  DESIGN 
L T2L1  A C3D  T L C3D 

NP 

T MASS  T 

LENGTH 

reactor 

i 

GAMMA  IC  GAM  EC  MIN' 

T1C2 

.020000 

.000969 

.068308 

4 

3.3061 

.12655 

.4923E  00 

.90R4C- 

06 

.4536E 

03; 

T1C3 

.020000 

.001139 

.078810 

6 

3.4973 

.13705 

.2095F  00 

.1024E 

07 

• 4686E 

os;. 

T1C4 

.020000 

.001308 

.088665 

6 

3.7102 

.14690 

• 19°7E  00 

.1051E 

07 

•4857E 

03 

T1C5 

.020000 

.001478 

.097868 

6 

3.9421 

.15611 

.1897E  00 

• 1 080E 

07 

,Fn48E 

n^ 

T1C6 

.020000 

.00164? 

.106433 

6 

4.1905 

.16467 

•1798E  00 

. 1 1 OAr 

07 

• 5255E 

T1C7 

.020000 

.001817 

.114374 

8 

4.4531 

.17261 

.9567E-01 

. 1 1 36E 

07 

.5470E 

07 

T1C8 

.020000 

.002262 

.132576 

8 

5.1065 

.19082 

• 8253E-8 1 

. 1 201E 

07 

.61 13E 

n 7 , 

T1C9 

.020000 

.002604 

.144293 

10 

5.8001 

.202*3 

.4718E-01 

. 1941E 

07 

.6644F 

03< 

T2C2 

.029200 

.000969 

.058624 

4 

4.0579 

.12606 

• 396°E  00 

• R752E 

06 

.4028E 

03 

T2C3 

.029200 

.001139 

.067846 

4 

4,2249 

.13520 

.3853E  00 

.9006E 

06 

.4169E 

07 

T2C4 

.029200 

.001306 

.076560 

6 

4.41 16 

. 1 4480 

• 1654E  00 

.9296E 

06 

.4374E 

03 

T2C5 

.029200 

.001478 

. 084769 

6 

4.61 60 

.15221 

• IS92E  00 

.9608E 

06 

.451 °E 

9\ 

T2C6 

.029200 

.001647 

.092480 

6 

4.8360 

.15QQ2 

. 1 52 RE  00 

.9936E 

06 

.4723E 

03 

T2C7 

.029200 

.001817 

.099712 

6 

5.0700 

.16715 

.1464E  00 

.1027E 

n7 

.4943E 

07, 

T2C8 

.029200 

.002262 

.116568 

8 

5.73R6 

.18401 

. 7290E- 0 1 

.1 1 14E 

07 

.5506E 

07' 

T2C9 

.029200 

.002604 

.127729 

8 

6.2971 

.19517 

.6631E-01 

. 1 1 74r 

07 

.61?9F 

0 3 

T3C2 

.040000 

.000Q69 

.052474 

4 

4.9889 

.13071 

.3405E  on 

.7952F 

06 

. 37R8E 

03- 

T3C3 

.040000 

.001139 

.060763 

4 

5.1318 

. 1 7000 

• 3321E  00 

.81  ROE 

06 

.3833E 

03 

T3C4 

.040000 

.001308 

.068613 

4 

5.3006 

,14685 

.3227C  nn 

• 8447F 

06 

. 7980E 

071 

T3C5 

. 040000 

.001478 

.076011 

6 

5.4955 

. 16425 

. 13RAC  00 

.P735E 

06 

.41 4RE 

03; 

T3C6 

.040000 

.001647 

. 082^78 

6 

5.6844 

. 161 2? 

• 1 340E  on 

.Q048C 

06 

.4334E 

03ii 

T3C7 

.040000 

.001817 

.089627 

6 

5 . 8Q63 

.16777 

.1291E  00 

.9S7GF 

06 

.45SfeE 

T3C8 

.040000 

.00226 2 

.184895 

6 

6.5837 

.1831 4 

.11 62E  00 

. 1026E 

07 

.FI 34E 

a 7 

T3C9 

.040000 

.002604 

.115152 

8 

7.0120 

.19339 

.6903E-01 

.1094E 

07 

.56F1E 

T4C2 

.050000 

.000969 

.048^74 

4 

5.8525 

.13721 

.30Q8E  00 

. 750  IF 

06 

.35^0E 

93; 

T4C3 

.050000 

.001139 

.056695 

4 

5.9Q4t 

.14493 

.302RE  00 

.7708F 

06 

.7ft44E 

03| 

T4C4 

.050000 

.001308 

.064004 

4 

6.1524 

.15224 

,?947E  00 

.794RF 

06 

.377RE 

03f 

T4C5 

.050000 

.001 476 

.070896 

4 

6.3256 

.15014 

.2R59E  00 

•8216E 

06 

.3931F 

03; 

T4C6 

.050000 

.001647 

.077369 

6 

6.5119 

.16561 

.122QE  00 

.8607E 

06 

.4102E 

os[ 

T4C7 

. OGOOOO 

.001817 

.083459 

6 

6.7102 

.17170 

.1187F  00 

.P813E 

06 

.U288E 

T4C8 

.050000 

.002262 

.097744 

6 

7.2782 

.18598 

• 1075E  On 

.9663F 

06 

.U842E 

03. 

T4C9 

.050000 

.002604 

.187305 

6 

7.7547 

.19554 

.992OE-01 

• 1 033F 

07 

.5325C 

03* 

03^ 

T5C2 

.060000 

.000969 

.046551 

4 

6.7324 

.14479 

.?R03E  00 

.7193E 

n6 

. 341 1 E 

T5C3 

.060000 

.001139 

.053863 

4 

6.8677 

.15210 

.2R2RE  00 

. 7,85E 

06 

.3516E 

03| 

T5C4 

.060000 

.001308 

.060779 

4 

7.0187 

.15902 

.2755E  00 

.7607F 

06 

• 3640E 

07!:. 

T5C5 

. 060000 

.001478 

.067293 

4 

7.1837 

.16553 

• 2674E  00 

.7854F 

06 

.3782E 

03 

T5C6 

.060000 

.001647 

.073414 

4 

7.3612 

.17165 

.259PE  00 

.01 23F 

06 

.7Q40F 

07^ 

T5C7 

.060000 

.001817 

.079148 

6 

7.5497 

.17739 

« 1 1 T 3E  00 

.0411F 

06 

.41 1 4E 

03g 

T5C8 

.060000 

.002262 

.0^2595 

6 

8.0889 

.19083 

• 1 0 1 1 E 00 

.°215E 

06 

.4631E 

T5C9 

.060000 

.002604 

.101576 

6 

8.5406 

.19982 

.9366E-01 

.9R55F 

06 

.5082E 

035; 

T6C2 

.060000 

.O0O96Q 

.043376 

4 

8,5852 

.1616? 

.2631E  00 

.6797F 

06 

,39*8E 

03^ 

T6C3 

.rroooo 

.001139 

.050130 

4 

8.6325 

.16837 

• 2572E  00 

. 6°67E 

06 

. 3353E 

93, 

T6C4 

.080000 

.001308 

.056505 

4 

8.7741 

.17474 

.2505E  00 

.7164F 

06 

• 3464E 

T6C5 

.080000 

.001478 

.062495 

4 

8.9285 

.18073 

.2433E  00 

.7304E 

06 

• 35°1 F 

93r; 
0 7 r ■ 

T6C6 

.080000 

.001647 

.0681 10 

4 

9, 0Q4t 

.18635 

.235RE  00 

•7622E 

06 

.3731E 

T6C7 

.OR0000 

.001817 

.073367 

4 

9.2698 

.19161 

.22R0F  00 

.7877F 

06 

.3885F 

03 

T6C8 

.080000 

.002262 

.085599 

6 

9.7700 

.20384 

.9234F-01 

.059RF 

06 

.4347E 

03 r 5 

T6C9 

.080000 

.002604 

.093725 

6 

10.1875 

.21196 

.0571E-O1 

.91 82E 

06 

.4749E 

03 

T7C2 

.100000 

.000R69 

.041342 

2 

10.2869 

.17958 

.9872F  00 

.6547F 

06 

.31 63E 

03P 

T7C3 

.inoooo 

.001139 

.047720 

4 

10.4001 

.18596 

.241 ir  00 

• 6705E 

06 

» 3252E 

03  . 

T7C4 

.100000 

.001308 

.083730 

4 

10.5852 

.19197 

•2347E  00 

.6886E 

06 

.?SC5E 

03';* 

T7C5 

.100000 

.001478 

.089367 

4 

10.6029 

.19761 

, 227RE  00 

• 7O07F 

06 

.3472E 

n7 ' 

T7C6 

.100000 

.001647 

.064637 

4 

10,8511 

.20288 

.22 07E  00 

.7^05E 

06 

• 3682E 

03  **v 

T7C7 

.100000 

.001817 

.069559 

4 

11.0185 

.20780 

•2135E  00 

.7537E 

06 

.3743E 

03  p 
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T7C8  .lnoooo  .002262  ,080958  6 1 

T7C9  .lOOOOO  .002604  .088481  6 1 

CERAMIC  AREA  ( C2  THRU  C9)S 
TAIL  LENGTH* 

total  ceramic  length  for  min  *zfoc*s 
total  MASS  (LBS.)  * 
total  lengths 
REACTOR  % 

PEAK  GAMMA  IC  $ 

MIN  • 7EOC * S 
END 


11.4937  .21020  .8645F-01  .81 96^  06  .4J65P 

11. 8888  .22672  .8028E-01  .R731E  06  .4573F 
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SECTION  3.3 

PLT  FAM  additional  sample  (linear-loq) 
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100 0.0  1250.0  1500.0  1750.0  5000.0  5250.0  5500.0  5750.0  6000.0 

FREQUENCY 


PLT  SRF  is  designed  for  3 dimensional  plotting  with  the  following  features: 

1)  Option  to  plot  a surface  (hidden  lines  are  suppressed),  or  families 
of  curves  in  3-D,  or  both. 

2)  Option  for  user  selected  Z axis  scales  or  automatic  Z axis  self 
seal ing  via  DISSPLA. 

3)  User  supplied  headings  as  well  as  X,  Y and  Z axis  labels. 

The  data  storage  scheme  for  X,  Y,  and  Z is  illustrated  in  Fig.  4.1. 

LPTS,  MEMBER,  and  NFAMLY  are  defined  as  before  where  now  each  member  corres- 
ponds to  a separate  value  along  the  Y axis.  Array  sizes  required  are  explained 
within  the  PLT  SRF  subroutine. 


\ not sm*  rwtro  -avra u/oc/90 


M/3W1  N.U.C.  OJJSPLrt  VEKSJO*  1 


i HorgQA  trwctto  tu'w w/ol/bo  aim a on  nr w sttt  >m  y 


SEADUCER  RUN  80122  HAY  6,1971 
IP  2.3  DESIGN  RUN  WITH  LOSSLESS  HEAD 
TRIAL  DESIGN  GENERATING  SCHENE  PO-5500 


MCCLFAPY.ftfttft^j; 


QRUN./RN  P|_TSRF  . 0 000  . LM60 1 07 1452?  . 2 » SO 

I3ASG.M  R,T,PLOT 
QASG.T  B..T.2004 
QREWIND  B. 

QCOPIN  P. . TPFS. 

QFRFE  R. 

GFOR.IS  CONTRL.  »CONTRL 

C*************  PLT  SRF  CONTROL  PmiTfMF  *********** 
COMMON  X^P(  fl),Y3D(  7).Z3D(  R.  7,  6),7MAT(  S6),W0PK3D(  T4)  , 

1 THDLPT  ( 6)  . IHDMEW(F,)»  TH0FAMM2,  6)  * THFAO( 1 ? . 3) » 

2 T BUF(2000)  c;rT 

COMMON  / MUMS  / LPTS*  MFMRFR.  NFAMLY,  MMLTmS 

□ IMP rjSION  I CHK(lR) 

DATA  M WRD  P c /l“/  D WORDS  PFP  C^Rf)  TMAGF  OM  IJMTVAC  '1AO 

DATA  IQ0LAR/6H  %/  110P 

DATA  IN  UNIT  / S / 


C**********  USER  RETS  N'HLlNS  = MUMRTR  OF  HEADTNR  LTMFS  TO  RF  LO^OFO 
NHL  I NS  = 3 

C**********  LOAD  IHFAD  FOR  PLOT  HEADINGS 
DO  4100  J = 1. NHL  INS 

RF  AO  ( IN  UMTT * 200  ) (IHFAD(T.J).  T=1,12) 

FLD  ( 30  » 6,  T HFAo  f 1.2  * J)  ) = I OOLAR  1 1 OP 

PRINT  ?0l  * ( IHF  AD ( T . J ) * T=1.12> 


4100  CONTimUF 


Qt  ********  * 

RFAD  AND 

print 

ANY 

EXTRA  INFORMATION 

TF 

• dfsirfh 

RFAD  ( 

IN 

UMTT. 

2o  0 ) 

(I 

CHK ( I ) . T = 1 » N WRD 

P 

C ) 

PRINT 

201  * 

(I 

CHK ( T ) . T = 1 . M WRD 

P 

C) 

RF  AP  ( 

IN 

UMTT. 

200  ) 

(I 

CHK (I) . T=1.M  WRn 

P 

C) 

PRINT 

201  . 

(I 

CHK (I) . 1=1. N WPD 

P 

0) 

200  FORMAT ( 13A6. A2)  1 1 OR 

201  format ( IX  * 13A6.a2)  iior 

C**********  IJSFR  LOADS  INPUT  DATA  FOR  LPTS.  ufmppr,  NF AML Y . X.  Y,  nTft  lRI 
POUNDS  = 2.20462 
LPTS  = R 
MFMRPR  = 7 
NFA^LY  - 6 
DO  R M = 1» MEMBER 
DO  P L = l * LPTS 

READ (5.100)  Y3D(m),  X3D (L).  (Z^n(L.M*N)  ,N=1 .MFflMLY) 

PRINT  101*  Y3D(m)  * X3D  (l_)  * ( Z3D  ( L * M * N ) * M=1  * NF  AMLY ) 

100  FORMAT  ( 3X  » 4X*3Ffl«G*3X*FR.4*FR.«sf  3E10.4) 

101  FORMAT ( IX  * 3X  * 4X*3FR.6*3X*FR.4*FR.Rf 3F10.4) 

DO  10  N = 1*NFAM|_Y 

IF  (M.EQ.2)  Z3D(l*m*N)  = 73D(L.M.N)  * POUNDS 
10  CONTINUE 

8 CONTINUE 

9 CONTINUE 


C**********  LOAD  IHD 
READ ( IN  UNIT* 
FLD ( 30  * 6.  IHD 
PRINT 

C***** *****  LOAD  IHD 
RFAD(  in  umtt* 
FLD < 30  * 6.  THD 
PRINT 


LPT  FOR  X AXIS  TITLF 

200  ) (IHD  LPT(I) . 1=1.6) 
LpT (6)  ) = I DOLAR 

201  • ( IHD  LPT(I) . T=1 *6) 

MEW  FOR  Y AXIS  TITLE 

200  ) (IHn  MEV(I) . 1=1 .6) 
Mem(6)  ) = I DOLAR 

20 1 » ( IHn  MEM( I) * 1=1 .6) 


1 1 OR 
1 1*8 


C**********  LOAD  THD  Fam  FOR  INDIVIDUAL  PLOTS  ( 7 AXIS  ) 
DO  4000  N = l.NFftMLY 
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I 


RF  AH  ( IM  I'MITr  ?f,  0 ) (IHH  PAM(I»M)» 

F L D < .10*  ft*  TMD  F A v { 1 P ► M ) ) = T DOL  AD  1 1 0 « 

PPIMT  2nl  * (IHH  F AM ( I » M ) r T = l,1?) 

4ono  continuf 

C * * * * * INPUT  IS  NOW  LOADED  * * * * * 

CALL  PLT  SRF 

c * * ♦ * * terminate  entire  plot  pun  ***** 

CALI  PLOT ( 0 . f 0.#  999) 

STOP 

END 


4-12 


DFOR,TS  PLTSRF,PLTSRF 

subroutine  PLT  srf 

C * * * * ♦ PLT  SRF  (USING  DTSSPLA  PACKAGE)  ************** 
C * * * * * GENERAL  PROGRAM  TO  PLOT  3D  c AML IES  OF  EURVFS  AND/OR  SURFACE'; 
C*******  PROGRAMMED  RY  L.  E.  MCCLEARY  NUC  CODE  601  ******* 


C 

c 

c 

c 

c 

c 

c 


* 

* 

* 


* 

* 

* 


1 

2 


* * * DIMENSION  AS  FOLLOWS  ***** 

* * * WHFRF  L=LPTS,  MrMEMREP.  N=NE AMLY , H=NHLTMS 

* LOAD  COMMONS  / /»  /NUMS/  FROM  CONTROL 

Common  X3D(L)*  Y3D(M),  730 ( L » M » N ) » 7MAT(L*m),  WORK  3D ( 2*L+P*M+4 ) 

IHDLPT(6)»  I HDMEM ( 6 ) * IHOFAM ( 1 2 » M) , THEAD(12»H) 
T BlJF(2f)0) 

COMMON  / MUMS  / LPTS,  MFMBER  * NF  AMLY  , NHLIMS 
Common  / MOLOAD  / 7MIN(N) » 7MaX(N),  ZPLOT(L),  YPLOT(L) 


i 

t 


COMMON 

1 

2 

COMMON  / 
common  / 


X 30 ( B ) » Y3D ( 7 ) , Z3D  ( «.  7>  6)*7MAT(  S6),W0RK3D<  34) » 
IHDLPT  ( 6 ) » IHDMEM  ( 6 ) » I HOF AM(1 2 » 6),  THEAD(12,3)r 
T 8UF ( 2000 ) 

NUMS  / LpTS , MEMBER,  NFaMLY,  MHLINS 
NOLOAO  / ZMIM ( 6 ) , ZMAX(6>,  ZPLOT(  B)  » YPLOT(  B) 


COMMON  / CM  VF3D  / XORIGfXTOP,  YORIG»YTOP,  ZORTG,ZTOP, 
1 XVU , YVU,  ZVIJ » XABS»YABS,ZABS 

DIMENSION  7 ORTG  L(  6),  ZOELTL ( 6) » ZTOPL ( 6) 

DATA  NBUF  / 2000/ 

DATA  IFTF/O/ 

DATA  I PLT  UN  / ft  / 


DATA (ZORTGL(L) 
DATA (ZD FLTL(L) 
DAT « (ZTOP  L ( L ) 


» Lr 1,6)/. 04 , 
,L=I ,6)/. 04, 
,L=l»6)/.20» 


S • , • 14  » 0 . » 

5.  , « 0 4 , .2, 

30 . , *26,  • B , 


.6F6,  300./ 

# ?Fft,  100./ 
1.4F6,  BOO./ 


print  1 
1 FORMAT (1H1) 

c * * * * * INITTALI7E  PLOT  BUFFER 

IF  ( IFTF  .EG).  0)  CALL  PLOTS  ( IRIIF , NPUF  , I PLT  UN) 
IFTrr 1 


c * * * * * FIND  max  Min  over  ALL  MEMBERS  for  facm  FAMILY 
□0  Soo  N = 1,  NFaMLY 
/ MTN(M)  = I.E30 
z MAx(N)  = -1.E30 
500  CONTINUE 

DO  3n00  L = 1,  LPTS 
DO  2000  M = 1,  IMfMBFR 
DO  1000  N = 1, NFaMLY 

if  ( Z3D (L ,M »M)  .LT.  Z MIN(N)  ) Z MTM(M)  = Z3B(L»M,N) 
IF  ( Z3D ( L , m , N ) . GT . Z MAX ( N ) ) Z MAX(M)  = Z3D(LrM,N) 

1000  CONTINUE 
2000  CONTINUE 
3000  continue 

c********  USE  nISSPLA  PACKAGE  ******** 

DO  5000  M = 1,  NFAMLY 

C**********  SFT  N SELF  7=0  FOR  USER  SUPPLIED  SCALES  ON  Z 

C**********  SET  M SELF  Z = 1 FOR  SELF  SCALING  ON  7 

N SELF  Z = 0 

N 5FLF  Z = 1 


* * 
* * 


SFT 

SFY 


SFT 

SFT 
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7 ORIG*  7 DPLT.  7 SCALP  ) 


r _ 

IP  ( M SELF  7 . Eq  , 0 ) GO  TO  20 
C * * * * * FTMD  7 SCAlPR  (OPTION) 

CALL  AXSPLT(  2 MlN(N) * Z MAX(N)*  a.* 
Z TOP  = 7 OPIG  + 7 DFLT  * 7 SCALP 

GO  To  30 


20  CONTinUF 

C**********  IJSPR  SETS  7 SCALPS  (OPTION) 

Z OPIG  = ZOPIGL(fj) 

7 oP|  T = ZOPLTL  ( N) 

7.  TOP  r 7J0PL(M) 

■^0  CONTINUE 

C*  ********  * USFP  SETS  X SCALPS 
X 0° IG  = .POOR 
X OFLT  = .noos 
X TOP  = .0028 

C**********  USFR  SPTS  Y SCALPS 

Y OPIG  = .0? 

Y DFLT  = .02 

Y TOP  = .10 

NT H PLT  = N 

CALL  RGNPL ( -MTH  PLT  ) 

CALL  TlTL^n(  JM  , t,  7.0*  7.0  ) 

do  4sno  j = i.mhltns 

CALt  HFADIM(  IMFaO(1»J)*  100.  2*  MHLINS) 

4500  CONTIMUP 
XARS  = 1. 

YARS  = 1. 

ZARS  = 1 . 

Call  AXFS30  ( IHn  L.PT.1O0*  THD  mEM.100*  IMP  PAM(  1 ,M)  * 1 on* 

A XARS,  YAPS*  ZAPS) 

X VI'  = -3. 

Y VI)  = -3. 

Z VI'  = 3. 

CALL  ViJARS(yVl).  YV'I*  7VU) 

CALL  GRAP3r(X0RIo*  XDFLT  * XTOP  * YOR T G * YOFLT , YTOP * 70RIG* 7DP!  T* 7T0P) 

DO  4qn0  M r 1,  yp'*p|ER 
DO  4300  L = 1*  1.PTS 

Y PLOT(L)  = Y30 ( M) 

ZPLOT(L)  = 73n(L,M»N) 

IP  ( M .MF.  1 ) GO  TO  4800 

C ♦ * * * * APPOw  PPRJrCTIONS  PDOM  CURVP  TO  AXTS 

CALL  PPI.PT.3  ( X’O  ( l ) * YPLOT  ( L ) * ZPLOT(L)*  XFROM,  VppOM ) 

CALL  RFLPT3 ( X3n ( | )*  YPLOT (L) » 70RIG*  X TO*  Y TO) 

CALL  QOAPOWf XPPOm.  YFROM,  X To*  Y TO) 

4800  CONTINUE 

C**********  IJSPR  SETS  N CURV  = 0 FOR  NOT  PLOTTING  ClIRVFS  OP  CONSTANT  Y 
C**********  USPP  SETS  M CURV  = 1 FOR  PLOTTING  CURVES  Oc  CONSTANT  Y 

N Cl'RV  = 0 
U CUpv  = 1 

IP  ( N CURV  .FO.  0 ) GO  TO  4820 
I MARK  = 1 

CALL  CURV3D ( X3D * YPLOT*  7PLOT*  LPTS*  T MARK) 

4820  CONTTMUF 
4900  CONTINUE 

CALL  VUPP30 

CALL  SURTPN( *POTh» ) 
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CALL  SUR5ZF(  Y3n(t)»  X3DILPTS)*  nr)(ll,  Y30 ( MFMBFD ) ) 

Lvi  r i 

00  4086  M = 1.  Mr^REP 
DO  4085  L = 1#  LpTS 
Z MAT(LM)  r Z3n ( i_ » M » M ) 

LM  = LM  + 1 
4985  CONTINUE 
49R6  CONTINUE 
IXPTOrl 
IYPT5=1 

CALL  5URMAT(Z  MAT*  IXPTS*  LPTS*  IYPTS*  MFMOE® . W0RK30I 
CALL  FMOPLt  -NTH  PlT  ) 

PRINT  49P0.  NTH  PLT*  ( THO  FAM(I*N)*  1=1*1?)* 

P YORlG*XnrLT,XTOP#  YORIG* YOELT  * YTOP*  70P TO * 7HFLT . 7T0P 
4990  FORMAT ( / • COMPLETEO  PLOT  M0.»,  13*  * ***  **  1?A6  / 
f * Y*,  3F 1 0 • 4 * * Y*»  3E10.4*  * 7* * 3^10.4  / ) 

5000  CONTINUE 

5000  continue 
rfturn 

END 


TFOP  » I S X3PMAT,y3nMAT 

FUNCTION  X3DMAT ( IX  ROW) 


COMMON  X3D( 


1 

2 


B)*Y30(  7)*?3r>(  B,  7,  B)»ZMAT( 

JHOLPT  f ft) # IHDMEm ( ft ) > iHnPAMf )?# 
T BUF ( 200 ) 
xsnnx  ROW) 


QFOR 


X3D«AT  = 

RETURN 

END 

IS  Y7dmat.y3dmat 

FUNCTION  Y7DMAT(yY  COL) 

COMMON  X3n(  B ) » Y30 ( 7)»730(  B » 7,  ft)#ZMAT( 

1 IHDLPT ( 6)  » lHDMEM(ft)»  THnPAM(12* 

2 I BUF (200) 

Y3DMAT  = Y3n(lY  COL) 

RFT"dn 

END 


BB ) t WOP^^D ( 34) 
ft>»  THEAD(12#3) 


Bfi),WORX3n(  ^4) 
f)»  THFAD ( 1 2 # 7 ) 


crT 


fTT 


CDF  OP,  IS  VIlPR^n,  V')Fr,3o 

SMRPOUTTMF  VUFR3r> 

COMMom  / Cm  VF^O  / XORIG,XTOP,  YO°TG» YTOP,  ZORTG.ZTOP, 
1 XVU,  YVll»ZVU»  X ARS  # YARS  » 7 ARF 


OIMFMSIOM  XOT(7),  YOT ( 7 ) » ZOT ( 2 ) , 
OTMFNSIOM  IX0TP(7»4),  TY0TC(7,4)» 
DATA  (ROFL(L) »L=1,3>  / -.003,  0., 
>ATA  ( ( IX0TC< I ,N) , 1=1 ,7) >N=1»4)  / 


YR(7),  YR(?)»  7R(?>»  POR  (3) 
I70TC(7  ),  lRXY7C(*»4) 

.003  / 


1,2, 7, 2, 1,1,1,  2, 2» 1,1, 1,2,1, 


OATA  ( ( IYOTC ( I ,N) , 1=1 ,7) »M=1.4)  / 
Y 2, 1, 1,1, 7, 2, 9,  1,1, 7, 2, 7, 1,2, 

OATA  ( IZOTCd  ) ,1  = 1,7)  / 

7 1 , l , 1 ,7,2,2, 1/ 

DATA ( ( IPXY7C ( I ,M) , 1=1 ,6) »N=1 ,4)  / 
R 1,7, 2, 3, 3,1,  7, 1,1, 3, 3, 7,  1,2, 


1,2,7,?,!  ,1,1  , 


2,1 ,1 ,1 ,2,7,2/ 
2,2,1 ,1 ,1,2,1/ 


1 , 

2,2, 3, 3, 1 

, 2,1 

,1,3, 3, 2,  1,2 

,2,3,3 

,1 

, 2,1, 1,3,3, 7/ 

IF 

( 

ZV'J  .LT, 

o.  ) 

PFT'JRN 

0 

COprjP  = o 

IF 

( 

XVU  .LF. 

0. 

.ANO.  yvu  ,LE • 

0.  ) 

M 

CORNP  = 1 

IF 

( 

XVU  .of. 

YAPS 

.AND.  yvu  .LE. 

0.  ) 

M 

rOPNR  = 7 

IF 

( 

XVU  .OF. 

X ApS 

.AMO.  YVU  ,GE. 

YAPS) 

M 

O0RMR  = 3 

IF 

( 

XVU  .Lr. 

0. 

.ANO.  YVU  .GE. 

YARS) 

N 

COPMP  r 4 

IF 

( 

N CORMP 

.fq.  o 

) return 

XOT ( 1 ) 
YOT ( 7) 
YOT ( 1 ) 
YOT(?) 
ZOT ( 1 ) 
Zn T(7) 


XHPIG 

XTnp 

yopig 

YTOP 

ZOPIG 

ZTPP 


00  700  I r'FL  = 1,3 

00  ?S0  I - 1 » G 

IC  = I XOT  C ( I ,NcORMP 1 

XR( 1 ) = YOT( IC) 

ir  = I YOTC ( I.MCOPNP) 

Y«(1  ) = YOT(IC) 

IC  = I 70TC ( I) 

Z«( 1)  = 70T  f IC) 

IF  ( I ,fq,  6 ) GO  TO  70S 
IC  = ixotc(t+i*ncormp) 

Xn ( ? ) = YOT(IC) 

IC  = lYOTC(!+l,MrORMP) 
yn(?)  = YOT(IC) 

I p = I70TC(T+1) 

ZP<2)  = ZOTfIC) 

50  To  210 
205  CONTy nUF 

IC  = I YOTC ( 1 , NfORNP ) 

X R ( ? ) = YOT(IC) 

IC  = I YOTC ( l.NrORMP) 

YR  ( 2 ) = YOT(IC) 
rc  = I70TC ( 1) 

ZP ( 2 ) = ZOT(IC) 

210  CONTTMIJF 

L WIDF  = IoyyZC(t,  MCOPNR) 

GO  To  (211,212,213),  L WT0F 

211  XR(1)  = XR(1)  + r OEL(T  OFL) 
XP ( 2 ) = XP( 1 ) 

GO  TO  PIS 

212  YR  ( 1 ) = YP(t)  ♦ n OEL(T  DEL) 


( ytop-xoptr) 

(YTOP-YOPTfi) 
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213 

215 

250 

260 


YR(?)  = YR(l) 

GO  TO  215 
ZR( 1)  = ZB(  1 ) 
ZR(?)  = 7B(1) 
CONTINUE 
CALL  CURV3D ( XP * 
CONTINUE 
CONTINUE 


+ r DEL (I  DEL)  ♦ (ZT0P-70RIG) 


YR»  ZB»  2 » 0) 


CALL  DASH 


IC 

- 

IXOTC  ( 

7»' 

YR 

(1  ) 

= XOT  ( 

IC) 

IC 

r 

IYOTC ( 

7 » 

YR 

(1) 

= YOT  ( 

IC) 

IC 

— 

IZOTC ( 

7) 

ZR 

(1  ) 

= ZOT( 

IC) 

DO 

300  I = t 

» 5 

IC 

r 

IXOTC ( 

I * 

XP 

(?) 

= XOT  ( 

IC) 

IC 

r 

IYOTC  ( 

I» 

YB 

(?) 

= YOT  ( 

IC) 

IC 

r 

IZOTC ( 

I) 

ZR 

(?) 

= ZOT  ( 

IC) 

CALL 

CURV3DI 

XBf 

NcORNR ) 
7 > NCORNR ) 


f 2 

I r NcORNR ) 
C) 

I » NCORNR ) 


Y«*  ZB*  2f  0) 


300 


CONTINUE 

call  RESET ( * DASH* ) 

RETURN 

END 
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QT  OP 

qprfp 

QPRT.T  TPFS. 

'aMAP.XI  SYMRtJM  » ARSRUN 

IN  CONTRL 

■axoT  absrum 

SFAQUCFR  RUN  oni 22  MAY  6»1971*  LIKE  80089  0=50 
IP  2.5  DESIGN  PUN  WlTw  LOSSLESS  HEAD* 

TRIAL  DESIGN  GENFRATInG  SCHEME  E0=5500  * 

CHARGE  TO  16000*01  ARRAY  DESIGN 


L T2L 1 

A C3D  T L C7Q 

NP 

T MASS  T 

LOjrTh 

ofactdr 

ga^ma  tf  gam  fc  mtm 

T1C2 

.020000 

.G00R6O  .088308 

4 

3.3061 

.126*5 

.49?3F 

on 

. Q984T 

06 

• 45*6F 

07 

T1C3 

.020000 

.001139  .078810 

6 

3.4973 

.13705 

.20O6F 

00 

. 1 o?4r 

07 

.46P6F 

07 

T1C4 

.020000 

.901708  .088685 

6 

3.7102 

,146°0 

• 1 997F 

no 

.1 051r 

07 

.4857F 

07 

T1C5 

.020000 

.00147R  .097868 

6 

3.9421 

.1*61 1 

• 1 897F 

on 

.1 naoF 

07 

.5048E 

O’ 

T1C6 

.020000 

.001647  .106433 

6 

4.1 Q05 

.16467 

.1708F 

00 

.1 108P 

n7 

• *?*5E 

07 

T1C7 

.020000 

.001817  .114374 

8 

4.4*31 

.17261 

.9567F-01 

. 1 1 36F 

07 

.5478E 

07 

T1C8 

.020000 

.002262  .132576 

8 

5.106* 

.1908? 

.B253F-01 

• 1 20 1 F 

07 

.61 13E 

O’ 

T1C9 

.020000 

.002604  .144293 

10 

5.8091 

.20253 

.471RF-01 

. 1241 r 

07 

.66t»4F 

O’ 

T2C2 

.020200 

.000969  .058624 

4 

4.0579 

.12606 

.3°69F 

on 

.8752F 

06 

. 4 028F 

07 

T2C3 

.02Q200 

.001139  .067846 

4 

4.2249 

.13529 

.3853r 

on 

.onofF 

06 

,t»169F 

07 

T2C4 

.020200 

.001304  .076560 

6 

4.4116 

.14400 

. 1654F 

on 

.92Q6P 

06 

. 4 ■*34F 

07 

T2C5 

.029200 

.001478  .084769 

6 

4.6160 

.15221 

.1*Q2F 

00 

. Q608F 

06 

. Ii5l  9F 

O’ 

T2C6 

.020200 

.001647  .092480 

6 

4.8360 

. 1*092 

. 1 *28F 

on 

. °°36P 

06 

.4723P 

03 

T2C7 

.029200 

.901817  .099712 

6 

5.0709 

.1671* 

. 1464F 

on 

. 1027F 

07 

.4943F 

07 

T2C8 

.020200 

.00226?  .116*68 

a 

5.7386 

.18801 

• 72QQF— 0 1 

.11 1«F 

07 

. C5«6F 

07 

T2C9 

.029200 

.002604  .12772° 

8 

6.2971 

.19*17 

.6631F-01 

. 1 1 74F 

07 

.61 29E 

07 

T3C2 

. 040000 

.O00°6n  .052474 

4 

4.980O 

.13071 

.3405F 

on 

• 7°5?F 

96 

. X708F 

07 

T3C3 

. 040000 

.901139  .060763 

4 

5.1318 

.13909 

. 3*21 F 

00 

.81  OOF 

96 

.3873P 

07 

T3C4 

.040000 

. O0130h  .068613 

4 

5.3006 

• 146*5 

.32?->F 

on 

. pU43F 

06 

. inflOF 

07 

T3C5 

.040000 

.901478  .076011 

6 

5.4853 

.1542* 

. 1 388R 

on 

• 8735F 

06 

. 4 1 U 8F 

03 

T3C6 

. O40000 

. HQ1647  ,08?Q78 

6 

5.6844 

.16122 

.134 OF 

on 

.00481- 

06 

. ti  774F 

07 

T3C7 

.040000 

.OOI8I7  .089527 

6 

5.8963 

.16777 

. 12°1F 

On 

,qX7PP 

06 

.U536F 

03 

T3GR 

.040000 

.002262  .104895 

6 

6.5037 

.18314 

.1 162F 

00 

. 1 026r 

07 

.51 3 4 P 

07 

T3C9 

.P40n0O 

.902604  .115152 

8 

7.012° 

.10339 

.6003F-01 

, 1094F 

07 

.56*1 F 

03 

T4C2 

.OR000O 

.000969  .048974 

4 

5.8*25 

.13721 

.3098F 

on 

. 750 1 P 

06 

• 3*30F 

07 

T4C  3 

.050000 

.901139  .056895 

4 

5.9941 

. 14403 

.3028F 

on 

.7708F 

06 

, "*644F 

03 

T4C4 

.050000 

.901708  .064004 

4 

6.1*24 

.15224 

,2°47F 

on 

• 7°48r 

06 

.3778P 

07 

T4C5 

.050000 

.901478  .070896 

4 

6.3256 

.15914 

.2859F 

on 

.8216F 

06 

.3°t1F 

03 

T4C6 

.050000 

.001647  .077369 

6 

6.51 1° 

.16561 

• 1 229F 

on 

• 8*07F 

06 

.41 0?P 

07 

T4C7 

.050000 

.901817  .083459 

6 

6.7102 

.17170 

, 1 1 87F 

no 

.8813F 

06 

.U288F 

07 

T4C8 

.050000 

.002262  .097744 

6 

7.2782 

.18*08 

. 1 07*F 

on 

.9663 P 

06 

• 4 84?F 

03 

T4C9 

.050000 

.902604  .107305 

6 

7.7547 

. 19**4 

.O929F-01 

• 1 933F 

07 

• 53?*F 

07 

T5C2 

.050000 

.90096Q  .046551 

4 

6.7324 

.14470 

.2893F 

no 

• 71 93P 

06 

,34  1 IF 

07 

T5C3 

. 080000 

.901139  .053863 

4 

6.8677 

.15210 

• ?828F 

on 

.7385P 

06 

.75I6E 

03 

T5C4 

.080000 

.001308  .06077° 

4 

7,0187 

.15902 

•?755F 

on 

.7607P 

06 

.3640F 

07 

T5C5 

.080000 

.001478  .067293 

4 

7.1837 

.16553 

.2674E 

00 

.7*54P 

06 

,778?F 

07 

T5C6 

.080000 

.001647  .073414 

4 

7.3612 

.17165 

.2590F 

on 

,81 23F 

06 

.■^rwinF 

07 

T5C7 

.060000 

.001817  .079148 

6 

7.5897 

.17739 

.11 13F 

on 

.8U1 ir 

06 

.4114F 

07 

T5C8 

.080000 

.002262  .092595 

6 

8.0«89 

.19083 

. 101 IF 

00 

• 92 1 *P 

06 

• 4 6 3 1 E 

07 

T5C9 

.080000 

.002604  .101576 

6 

8.5406 

.19982 

.9366F-01 

, 98*5p 

06 

.*nP2F 

o7 

T6C2 

.080000 

.000969  .043376 

4 

8.5052 

.16162 

• 2631 F 

no 

• 6’?97F 

06 

.32*8E 

07 

T6C3 

.080000 

.001139  .050130 

4 

8.6325 

.16837 

.2572F 

on 

. 6°67P 

06 

. 3353F 

07 

T6C4 

.080000 

.001308  .056505 

4 

8.7741 

.17474 

. 2*05F 

on 

.7164P 

06 

. 3464F 

07 

T6C5 

.080000 

.001478  .062495 

4 

8.9285 

.18073 

.2433F 

on 

.7384F 

06 

.3501F 

03 

T6C6 

.080000 

.901647  .068110 

4 

9.0Q41 

.18635 

.235BF 

00 

.76?2P 

06 

.3731F 

03 

T6C7 

.080000 

.001817  .073367 

4 

9.2698 

.1°161 

.2280F 

on 

.7877F 

06 

. T885F 

07 

T6C8 

.080000 

.002262  .085599 

6 

9.7700 

.203*4 

.9234F-01 

#ppoor 

06 

.4347E 

07 

T6C9 

.080000 

.002604  .093725 

6 

10.1875 

.21196 

.8571F-01 

.91 8?p 

06 

.4749F 

0’ 

T7C2 

. 1 00000 

.000969  .041342 

2 

10.286° 

.179*8 

.°87?f 

on 

.6*47F 

06 

.31 63F 

07 

T7C3 

.100000 

.001139  .047720 

4 

10.4093 

.18*06 

. 241 1 F 

00 

• 6705F 

06 

• 3752E 

07 
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T704 

. 1 oooon 

.O013Or 

.053730 

4 

10.545? 

.19107 

.2347E 

00 

.6886F 

06 

.7355F 

07 

T7C5 

.100000 

.00147ft 

. 059367 

4 

10.6^29 

.19761 

.227RE 

00 

.7087F 

06 

.3472F 

07 

T7C6 

. 1 00000 

.001647 

.064637 

4 

10.8511 

.?02ftfl 

.2207F 

00 

• 7305F 

06 

.3602F 

03 

T7C7 

.100000 

.001917 

.069559 

4 

11.0185 

,?0780 

•2135F 

00 

.7537F 

06 

.3743F 

03 

T7C9 

.1 00000 

.002262 

.080959 

6 

11.4937 

.21920 

.8645F-01 

.8196F 

06 

.4165F 

03 

T7C9 

.100000 

.002604 

.088*81 

6 

11.8«85 

.22672 

.ft02flr-01 

.8731F 

06 

.4533F 

03 

CERAMIC 

area 

( C2  THRU 

C9>$ 

TAIL  LENGTH  t 

TOTAL  CERAMIC  LENGTH  FOR  MIN  »ZEOC»* 
TOTAL  MASS  (LBS.)  -5 
TOTAL  LENGTHS 
REACTOR  s 
PEAK  GA^mA  IC  * 

MIN  *ZFOC»  1! 

QFIN 
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SECTION  4.3 

PLT  SRF  additional  sample 


4-21 


SECTION  5. 
PLT  PTN 


PLT  PTN  is  a modified  version  of  PLT  FAM  specifically  oriented  towards 
polar  plot  applications  (such  as  beam  patterns).  Its  features  include: 

1)  Option  for  user  selected  axis  scales  or  automatic  self  scaling  via 

DISSPLA 

2)  User  supplied  headings  and  legend  labels. 

The  data  storage  scheme  for  angle  vs.  magnitude  is  the  same  as  for 
Fig.  3.1  of  PLT  FAM  where  X corresponds  to  angle  and  Y corresponds  to  magnitude. 
Array  sizes  required  are  explained  within  the  PLT  PTN  subroutine. 

Notice  that  in  this  example  the  angular  orientation  has  been  rotated  90° 
so  that  0°  is  verti cal . 


SECTION  5.1 

PLT  PTN  plotting  sample 


5-2 


90  SOU  3190  I OW  10TJ  80KI  (X  8J 


SECTION  5.2 
PLT  PTN  program  listi 
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2000 


601 

M3S=DL0T 


M,nnn*n»  no.  03 


L MO0LrARY  P|  T otw  hcam  pat-' 


jor 

MOUNT 
PORTP Am 

***********  PLT  pTK'  CONTPOL  PfMJTT»jF  *********** 
COmyom  / / XDAT A ( 361 ) » YDATA(361»  4,1) 

COMMON  / CMME AD  / T HIDMEM  ( 6)  , T Hr  AO  (12,3) 

COMMON  / fH  'MS  / LPTS,  Mfmrfr,  MFAMLY , MMLTNS 
COY  mon  / CPI.  LPL  / DTA  LRL  ( 4,2),  IPK  l Rt.  ( 4, 12) 

U TMFvjsI ON  T MOPE  ( 3 ) 

DATA  IDOL  AP/5H  % / P^n 

DATA  IM  UNIT  / 5 / 


C**********  IJSPp  SETS  NHLINS  = Mliyppp  OP  HEADTM6  LTA'FP  TO  RF  i.oahfp 
NHL T MS  = 3 

C**********  LOAD  THFAD  pop  plot  HFAPIMGS 
DO  UioO  J = 1 » NHL T NS 

RFAD(  TM  UNTT,  2(10  ) (IHFAD(T,J),  T=1,12) 

CALL  SFLD(  24#  6,  T HFAD(12,J)»  TnoLAP  ) P’O 

PPIMj  201  r ( IHr  AO ( I , J ) , T — 1 ,12) 

umo  COMTTNUF 

200  FORMaT(16AS)  1230 

201  FORMAT  ( 1 y r 1 GAS)  I^P 

r**********  | OAn  THn  mrm  FOP  LFOFND  TTTl.F 

C'rAn(  im  UNIT,  200  ) f I HD  MFm(T),  Trl,6) 

CALL  SFLO(?4,  G , ThD  MFM(a),  IDOLAR)  1^30 

PP I MT  201  , ( I HD  MFM(I),  T- 1 , 6 ) 

N TImfS  = 1 

mp  S000  I T T MFS  = 1»  M TIMES 


r**********  HSFR  LOADS  INPUT  DATA  fop  LPTS  , MFmdcd,  MFAMLYf  Yf  V,  ATf  ( n, 

ifmpfr  - 1 
MFA"LY  = 1 

lirflA(Si°Q)  MAMO»ST  AMGfDAMR,  e-RFr)  , THAMG  » TMODF  (1  ) , TMODF  < 2 ) • TMpnr  ( 3 ) 
PRINT  ORf  MAMR  f ST  ANGfDAMGf  frfo  , IHAMP  , IMODF  ( 1 ) , TMOnr  ( ?.)  r Tmoof  ) 
oq  FOPMftT ( IlOf  3F 1 0 . 2 f SX.A5,  SX,3AS) 

LPTS  = NANG 
m r i 
J = 1 

JTA  LRL (Mr  1 ) = FpFO 

DO  0 l - If  LPTS,  10 
LLO  = L 
LHI  = L+P 

IF  ( lmt  .GT.  NAmG  ) LM I = NANG 

iIF AO ( S » 1 0 0 ) (Y  DATA(LLfM,M) , LL=LI  0 , LMT ) 

PRINT  100,  (Y  0 aT  A ( LL » M » M ) , LL=LL0 , L^T  1 
100  FDRM/vTdOFR.l  ) 
n COajt tnUF 

R QUAD  = on. 

X = ST  AMG  + P QIJAD 
D01P  L = 1 f LPTS 

I F ( y DATA (L, M»N)  .LT.  -3P,5  ) Y 0ATA(L»m,N)  = -*o.S 
X DATA(L)  = X * -A.1415P  / 1R0. 

X = Y + HAND 
18  comttmuf 

C * * * * * INPUT  IS  NOW  I.OAOED  ***** 

CALL  PL.T  PTM 


U 


SO no  COMTTNUF 
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S"RDOtlTIMF  PLT  PTN 

C * * * * * PLT  PTfj  (USING  D TSSDL A PACKA<5r)  ************41*** 

C ♦ * * * * GFNERAL  PROGRAM  TO  PLOT  FA^TLTFR  PF  CURVE?  WTTM  LFGFNDS  * * * * 

C*******  PROGRAMMED  PY  L.  Ft  MCCLEAPY  Ml  IF  CODE  ^,0 1 ******* 

C * * * * ♦ DIMENSION  A?  FOLLOWS  ***** 

C * * * LO AO  COMMONS  / /»  /CMHFAP/,  /Nl|MS/»  AMD  ARRAY  nT A I.BL  IM  CONTROL 

C * ♦ * * * WHFRF  LrLPT?*  m=m£moFR»  N=NEAMLY,  H=MHL T NS 

C COMMON  / / XDATA(L)*  YOATA (L»M»M) 

C COMMON  / CMHFAO  / TH0MFm(6)»  IHFAO<1?»m) 

C COMMON  / MUMS  / LPTS»  MFMqpo,  MEA«LY»  ABLINS 

C COMMON  / CPL  LRL  / OTA  lRI.  ( M,?),  IPK  L°L  ( M*1S>> 

C COMMON  / NO  LOAD  / YMIM(N)*  YmAX(N)»  YDLOT  (L) 

COMMON  / / XOATA(3M)»  YOAT A ( 361  » tt»l) 

COMMON  / CMHFAO  / THDMCM(  ft)  » lHrAD(19,3) 

COMMON  / MUMS  / Ldt?»  Mrwpfp,  MR  AMLY  t MHLTMS 
Common  / CRL  LRL  / OTA  LRL ( 4,?),  IPK  LRL  ( If  12) 

COMMON  / NO  L0«0  / YMJN  ( 1 ) t Y MAX  ( 1 ) » Y P|_0T(3ftl) 

DTMFNSION  Im  LRL  m ( 4 ) 

COMMON  / Cm  PLOT  / I Rl  IF  ( 200) 

DATA  MBUF  / 200  / 

DATA  iFTR/n/ 

DATA  I PLT  UN  / ? / 

C * * * * + INITIALI7F  PLOT  PUFFER 

IF(IFTF  .FO.  0)  CALL  PLOTS  ( I RUF  t NRIJF  t T PLT  UN) 

IFTFrl 

C * * * * * FINO  max  MIN  OVER  ALL  MEMBERS  FOP  EACH  FAMTLV 
DO  S00  N = It  NFaMLY 

Y MIN(N)  = l.F^O 

Y MAX(N)  = -l.F3o 
500  CONTINUE 

DO  3000  L = 1 # LPTS 

no  ?noo  Mr  it  member 

DO  1000  N = 1 t MFaMLY 

I R ( Y DAT  A f L » M t n ) .LT.  Y MTN(N)  ) Y MTN(M)  = Y D6TA(I_*M»M) 

IF  ( Y DATA(L»m^|)  ,OT.  Y maX(N)  ) Y maX(n)  = Y r)ATA(L*M»w> 

)000  continue 
2000  continue 
3000  continue 

c ********  USE  DTSSPLA  PACKAGE  ******** 

DO  FnOO  M = 1»  NFAMLY 
CALL  RGNPL.(  N ) 

M = 1 

DO  4400  L = 1 t LPT? 

4400  Y PlOT(L)  = Y OATA(L»M»N) 

C**********  USFR  SETS  I SCALP  = 0 rOR  SELF  ?CAl  TM© 

C**********  USFR  RETS  I SCALE  = 1 FOR  USFR  SURPl.lEn  ?CAL.FS 

I SCALE  = 0 
I SCALE  = 1 

I F ( I SCALP  .EO.  0 ) GO  TO  4440 


«;rT 

crT 

err 


c**********  FOR  I'SEP  SUPPLIED  scaler t 
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£**********  tJSFR  SETS  R ORIGr  R STEP  * R INCH  FOR  PACH  FAMILY 
R OPIG  = -40, 

R 5TFP  = 10. 

R INCH  = 4, 

30  To  4445 

4440  CONTINUE 
C * * * * * SELF  SCALING 

CALL  AXSPLT(Y  miN(n)»  Y MAX(N)*  4.*  R PRIG,  p STFP,  R TNCH) 

4445  CONTINUE 

C * * * * * SCALFS  ARE  NOW  AVAILABLF  FOP  PLOTTING 
R TOP  = P ORIG  + R STEP  * R INCH 
CALL  P RANGE ( R ToP) 

CALL  POLORG ( R ORlG) 

0 INCH  = P INCH  + R INCH 
ITITLE  = 1 

CALL  TITLE ( 1H  *ItTTLE*  1H  *1,  1H  *0,  D INCH»  D INCH) 

TH  CONV  = 1, 

CALL  POLAP(TH  COnV,  R STEP*  R INCH,  R TNCH) 

00  4S00  J = 1 • NHL I NS 

4500  CALL  HEAOIN ( IHEaD(1*J)*  100*  2*  MHLINS+1 ) 

CALL  HEADIM(1H  *1*  2*  MHLTNS+1  ) 

CALL  GRIO ( 1*1) 

CALL  MARKER(  V ) 

1 MARK  = 0 
I MARK  = 1 

CALL  CURVE ( X OATA*  Y PLOT*  LpTS*  I MARK  ) 

C**********  USER  SELFCTS  LEGEND  FORMAT  FOP  ENCODF 
4 FORMAT ( F7 . 1 , • ** ) 

ENCODE  1 20  * 4*  TH  LBL  M)  DT  ALBL (M  * 1 ) 

CALL  LINES ( TH  LBL  M*  IPK  LBL*  M ) 

IP  ( mFMRFP  ,E0.  1 ) GO  TO  4901 
DO  4900  M = 2 * MEMBER 
00  4600  L = 1 * LPTS 
4600  Y PLOT ( L ) = Y DATA(L*M*N) 

CALL  MARKPR(M) 

CALL  CURVE(  X DATA*  Y PLOT*  LPTS*  I MARK  ) 

ENCOnF ( 20  * 4*  TH  LBL  M)  DT  ALBL (M  * 1 ) 

CALL  LINES ( TH  LBL  M*  IPK  LBL*  M ) 

4900  CONTINUE 

4901  CONTINUE 

C**********  USER  SETS  N LEGNO  = 0 FOR  MO  LEGEND 
C**********  USER  SETS  N LEGND  = 1 FOR  LFGEND 
N LFGND  = 0 
N LEGND  = 1 

IE  ( N LEGNn  ,EO.  0 1 GO  TO  4960 
XLEGND  = D TNCH  - 1.5 
YLEGND  = 0 TNCH  + .5 

CALL  LEGEND  (IPK  LBL*  MEMBFR , XLEGND*  YI.EGND  ) 

CALL  MESSAG ( IHD  M£M*  100*  XLEGND*  YLEGN0-.2  ) 

4960  CONTINUE 


CALI.  HFTKHTf.in) 

Call  ^FSSAG(tn  nc-fiw  «;»  p T''Ch*.s-»  i n»  n iMru+.joi 
CALL  AMGLF(oo.) 

CALL  MFSFAGf *oo  nFr,',  6,  -,it  n TMCH*.r-.9«;i 
CALL  pFSFTI • AMGLF* ) 

CALL  RFSFT( »MFTGmT* ) 

CALL  FNDPLf  tl  ) 

S'lno  COMT  TMtIF 


RFT"rm 

INC 
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f. 


■SEND 

SDATA 

MARCH  6.1971* 

BEAM  PATTFRN  THFTA  * 
TR ANSPARFnT  * 
FREQUENCY* 


46 

Rn.on 

1 

no 

• 

o 

o 

5^0.00 

THFT  A 

TRANSPARENT 

.n 

-.0 

-.1 

-.3 

— • 4 

-1.2 

-1  .7 

-2.7 

-7.- 

-3.4 

-4.1 

-4.9 

-5.4 

-6.4 

-7.4 

-9.0 

-lo.-* 

-1  1 .4 

-17.1 

-IS. 2 

-17.7 

-19.7 

-22.4 

-24.9 

-26.9 

-25. 0 

-23.7 

-71 .4 

-79.  • 

-19.  n 

-16.1 

-17.3 

-16.7 

-16.7 

—15,5 

-16.4 

-16.1 

-14.9 

-14.  t 

-14.  S 

-14.3 

-14.1 

-14.0 

-1  3.9 

-13.  P 

.0 

.9 

.9 

PND 


I 

I 


L 


jt  N) 


SECTION  5.3 

PLT  PTN  additional  samnle 
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